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Abstract 

The WMAP results on the scalar spectral index n and its running with scale, 
though preliminary, open a very interesting window to physics at very high en- 
ergies. We address the problem of finding inflaton potentials well motivated 
by particle physics which can accomodate WMAP data. We make a model 
independent analysis of a large class of models: those with flat tree-level po- 
tential lifted by radiative corrections, which cause the slow rolling of the infla- 
ton and the running of n. This includes typical hybrid inflation models. In 
the small-coupling regime the predictions for the size and running of n are re- 
markably neat, e.g. —dn/d\nk = (n — l) 2 <C 1, and n does not cross n = 1, 
contrary to WMAP indications. On the other hand, n can run significantly if 
the couplings are stronger but at the price of having a small number of e- folds, 
N e . We also examine the effect of mass thresholds crossed during inflation. 
Finally, we show that the presence of non-renormalizable operators for the infla- 
ton, suppressed by a mass scale above the inflationary range, is able to give both 
dn/d\nk ~ O(-0.05) and N e ~ 50. 
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1 Introduction 



The measurement and analysis of the cosmic microwave background (CMB) by the 
WMAP collaboration [1-3] has provided results of unprecedented accuracy on the 
spectrum of the CMB anisotropies. The data are basically consistent with the main 
features of inflationary cosmology [3], implying important constraints for particular 
models (many of which can be discarded). One of the most intriguing WMAP results, 
though still uncertain, indicates that the spectral index of scalar perturbations, n, may 
exhibit a significant running with the wave-number, k. More precisely, assuming a 
linear dependence of n with Ink, n{k) = n(k ) + [dn/dlnk] ln(/c//c ) {k Q is some chosen 
"pivotal scale", k = 0.002 Mpc -1 ), the analysis of ref. [3] gives n(k ) = 1.13±0.08 and 
dn/d\nk = -0.055^;^ at 68% C.L. (These values are obtained from WMAP data 
combined with different large-scale-structure data.) This means in particular that 
n — 1 passes from positive to negative values in the range of ^-values corresponding 
to these data (10~ 4 < k Mpc < 4). It should be said, however, that a fit assuming 
constant spectral index is also consistent with the data, although somewhat worse, 
giving n = 0.99 ± 0.04 (with WMAP data only [1]). 

As a matter of fact, most inflationary models [4] do predict that n runs with k, but 
the details of the "observed" running are very difficult to fit [5]. The last scales entering 
the observable universe (which we will denote k* ~ 10~ 4 Mpc -1 ) should have crossed 
the horizon ~ 50 — 60 e-folds before the end of inflation, at a time we will call t*. The 
range of k probed by the WMAP data corresponds to the first ~ 6.5 e-folds after 
(this range extends to ~ 10.5 e-folds if one adds large-scale-structure data). During 
this "initial" period of inflation, n — 1 should change sign (if the WMAP indication is 
confirmed), and it is not easy to cook up an inflationary potential implementing this 
feature (see [5,6] for an incomplete list of several analyses along these lines). In the 
present paper we show, however, that a large and very interesting class of inflationary 
potentials, including those of many supersymmetric hybrid inflation models, are in 
principle able to accommodate this behavior quite naturally. 

Let us recall a few basic results of inflationary cosmology in the slow-roll approxi- 
mation. The slow-roll parameters are given by 
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where 8ttM 2 = 1/GNcwton, V is the scalar potential of the (slow-rolling) inflaton field, 

0, and the primes denote derivatives with respect to <fi. The slow-roll approximation 

requires e, \q\, |£| <C 1, and the failure of this condition marks the end of inflation. The 

number of e-folds produced during the slow-rolling is given by 

/■tend 1 r<P V 1 r<t> 1 

N e = Hdt~— —M = — / -^d<P , (2) 

Jt Ml J<p cnd y y M p J<p cnd ^ ^ K J 

where H is the Hubble parameter [with H 2 = V/(3M 2 )], meaning that the scale 
parameter of the universe, a, grows as a — > e^a. Successful inflation requires 1 N e > 
50 — 60. At first order in the slow-roll parameters [8], the power spectrum of scalar 
perturbations, Pk (called An(k) = (2.95 x 10~ 9 )A(k) in ref. [3]), and the spectral index, 
n, are given by 

1 V 
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(3) 



n,l + ^l + 2,-6e. (4) 

Let us remind that, in the slow-roll approximation, the derivatives with respect to In A; 
can be related to field derivatives by using 

as is clear from (2) and aocl/k. It is then trivial to obtain 

d " ~ -2£ + 16e V - 24e 2 . (6) 



d\nk 

The WMAP analysis [2,3] gives 

P k = (2.95 x 10~ 9 ) x (OJOljui) at k = k = 0.002 Mpc' 1 , (7) 

and n as mentioned in the first paragraph: 

n(k ) = 1.13 ± 0.08 , ^ = -0.0551°;^ • (8) 

Since usually e 77, then n ~ 1 + 2i], and the change of sign in n — 1 must arise from 

a change of sign in rj (which implies in turn that V must pass through an extremal 

1 It has been argued, however, that the required 50 — 60 e-folds could be achieved through several 
steps of inflation, corresponding to different inflaton fields and/or inflaton potentials [7]. WMAP 
would only be sensitive to the first (s) of these steps. 
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point), something which, as mentioned above, is non-trivial to implement (see ref. [5] 
for a systematic analysis of this and related issues). 

On the other hand, the evidence for a running n comes from the first (I < 5) 
WMAP multipoles [9], which corresponds to the first ~ 1.5 e-folds of inflation (after 
t*). Excluding those points, the fit is consistent with a constant index, n ~ 0.94±0.05, 
i.e. the most significant running of n seems to be associated to the first few e-folds 
and then n remains quite stable. In any case, notice that the errors associated to these 
analyses are big (actually refs. [3] and [9] are only consistent within these errors). 

Hence the results concerning the spectral index can be summarized by saying that 
they indicate 

3E* = o <- -<*>. < 9 > 

during the first few e-folds of inflation after t*, and 

n~ 0.94 ±0.05, (10) 

during the following e-folds. Of course, if the previous dramatic running is not con- 
firmed by forthcoming analyses, n could still exhibit a milder running, compatible with 
(10) in the whole range of observed scales. 

In section 2 we study a large class of inflationary models defined by having flat 
tree-level potential, with radiative corrections being responsible for the slow-roll of 
the inflaton and the running of n. The results depend on whether the couplings are 
weak (subsection 2.1) or strong (subsection 2.2) but generically one cannot get both a 
significant running of n and a large N e . We devote section 3 to a numerical analysis of 
this issue in the well motivated type of D-term inflation models. Things improve in this 
respect if inflation is affected by some high energy threshold, possibility which is studied 
in section 4. Subsection 4.1 studies what happens if the inflaton crosses such high- 
energy threshold shortly after t* and subsection 4.2 analyzes the most promising case 
in which the threshold is above the inflationary range but affects inflation through non- 
renormalizable operators. In this last case we can easily accomodate both dn/dliak ~ 
0{— 0.05) and N e ~ 50 — 60. Finally, Appendix A gives some technical details relevant 
to section 3 and subsection 4.1. 
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2 Spectral index in a large class of inflation models 



Here we analyze a large class of inflationary models, namely those with a flat inflaton 
scalar potential at tree-level. Then, the flatness is spoiled by radiative corrections, 
which are responsible for the slow-roll of the inflaton and, in turn, for the value and 
running of the spectral index. This class of models include many supersymmetric 
hybrid inflation models as a particularly interesting subclass (for other types see e.g. 
[4,10]). To be concrete, we will focus in the following on hybrid inflation models, but 
all the results hold also for generic tree-level-flat inflation models. 

Hybrid inflation models [11] are among the most popular and interesting ways to 
implement inflation. Their main characteristic is that there are two (or more) scalar 
fields involved in the inflationary process: the slow-roll field, 0, and the "waterfall" 
field(s), H. The scalar potential is such that, as long as is beyond a certain value, C , 
the minimum of the //-potential is at some "false" vacuum (typically at H — 0) and 
has an exactly flat potential at tree-level. Below C the minimum of the //-potential 
is at the true vacuum. Once H departs from the origin, the potential is not flat 
anymore, inflation ends and goes to its true minimum as well. The slow-roll of is 
driven by the radiative corrections to the potential, which lift the tree-level flatness. 
The main merit of such hybrid inflation is that the slow-roll parameters e and r\ can 
be easily small, as required, due to the loop-suppression factors and the fact that the 
scale of is not directly related to the magnitude of V. 

For our purposes, the main feature of hybrid inflation is that during the inflationary 
period the tree-level potential, V = p = constant, is exactly flat for 0, all its curvature 
arising from radiative corrections. Generically, at one-loop 

V(0)=P + /? ln^, (11) 

where Q is the renormalization scale (which might have absorbed finite pieces). Note 
that p depends implicitly on Q through its renormalization group equation (RGE) and 
that the Q-invariance of the effective potential implies 

at one-loop. Finally, m(0) is the most relevant 0-dependent mass [there could be several 
relevant (and different) masses, which is a complication we ignore for the moment]. 



4 



The leading-log approximation (which amounts to summing the leading-log contri- 
butions to all loops) is implemented in this context by simply taking Q = m(<f>), so 
that 

V(<j>) ~ p(Q = m(0)) . (13) 

In general, one expects m 2 (0) = M 2 + c 2 2 , where M does not depend on and c 
is some coupling constant. We will ignore for the moment the possible presence of 
the M piece, as indeed happens in typical hybrid inflation models. From the equation 
Q = m(4>) we see that each value of Q corresponds to a particular value of through 2 
= Q/c(Q) (note that c also depends on Q according to its own RGE). 

The 0-derivatives of the potential (13) and the corresponding slow-roll parameters 
can be easily calculated using 



dQ Q 

— — = — a = cot 

d(p 



where 



In particular, 



a c 



and P c 



dc 



d\nQ 



(14) 



(15) 



V = 
V" = 



(3 





— a 
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(16) 



where the dots stand for derivatives with respect to In Q. 

Notice that from eq. (5) we can relate the wave-number k with the scale Q: 

3H 2 



, k 
In — ~ 

k 



P 



-XQi-Q' 



(17) 



(3c 2 M 2 ^ ' (3c 2 

(This expression is exact in the small-coupling approximation.) Therefore, the spectral 
index at different k scans physics at the corresponding (high-energy) scale. We analyze 
next this connection in detail. 



2.1 Small-coupling regime 

In the regime of very small coupling constants one has (3 c /c 1, and thus a ~ 1, and 
we can also neglect all terms within the square brackets in (16), except the first one. 
2 We take c and <j) as positive, which is always possible through field redefinitions. 
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Then 

so e C ?], as usual, and thus n — 1 ~ 2r/. 

From the previous expressions it is clear that in this regime there exist very strong 
and model-independent relations between the observable parameters. The running of 
n with In k is given by 

5^ = -2£»-V —(«-!)'. 0») 

where we have used eqs. (6) and (18). Consequently, dn/dhik is negative, as suggested 
by observation, though its value tends to be quite small. In fact the sign of n — 1 cannot 
change along the inflationary process, which would contradict the WMAP indication 
of n crossing the n — 1 value. The sign of dn/d\nk and the correlation of its absolute 
value with that of n — 1 are quite model-independent: they apply whenever the field 
dependence of the inflaton potential is well described simply by a In 0. If the value of (3 
changes along the inflationary course, e.g. if some threshold of extra particles is crossed 
in the way, this correlation still holds except in the neighborhood of the thresholds. 
The signs of the first derivative of the potential, V ~ /3/<f>, and 77 (and thus n — 1) are 
also correlated; namely 

sign(77, — 1) = — sign(y') = — sign(/3) . (20) 

Notice that a change in the sign of n — 1 can only occur if V vanishes, indicating the 
presence of a local minimum. Usually (3 (and thus V) is positive and therefore we 
naturally expect n < 1. Then rolls towards the origin, where the true minimum of 
the (b— potential is in typical hybrid inflation models [11]. Inflation ends either because 
the flatness condition gets violated below some critical value C (as mentioned above) 
or because e or (much more often) \q\ become 0(1). In this respect, notice that e and 
1 77 1 increase as rolls to smaller values. On the contrary, if (3 (and thus V) is negative, 
then n > 1, runs away and the slow roll parameters e and 77 get smaller as time goes 
by. Inflation would not stop unless reaches some critical value where the tree-level 
flatness ends and all the fields go to the true minimum. This is not the most common 
situation, but it cannot be excluded (see e.g. [12-14]). 
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In the usual case, with inflation ending because 77 gets 0(1) (which according to 
the previous discussion requires rj < 0) and in the absence of thresholds of new physics 
during the inflationary process 3 the number of e-folds, N e , since t* until the end of 
inflation can be easily computed plugging (18) in eq. (2) and is simply 



N e ~ -- 
2 



2??(0*) 



21) 



If there are no other steps of inflation (see footnote 1), then N e = 50 — 60; otherwise 
it could be smaller. Then the value of the spectral index at the beginning of inflation 
is given by 

ra„ = n{K) ~ 1 + 277(0,) - 1 " , (22) 

which means n* < 0.983 (the equality holds for N e = 60, i.e. if there are no additional 
steps of inflation). For a given value of (3/p this also fixes the initial value 0* through 
eq. (18). In addition, we can use (22) to write eq. (19) in an integrated form as 

H = 1 " N e - ln(k/K) ■ (23) 

For later use, it is convenient to rewrite this in terms of the initial value of the spectral 
index, n* given by eq. (22), as 

11 

n ~ U * + W e ~ N e -\n(k/h) ■ { > 

If inflation ends by the breaking of the flatness condition, eqs. (22-24) hold, replac- 
ing N e — > — 1/[2?7(0*)]. This also holds for the (3 < (and thus n > 1) case. Finally, 
note that there is no problem in reproducing the observed scalar power spectrum (7): 
for a given value of fi/p, this can be achieved with an appropriate value of V = p, 
according to eq. (3). 

In summary, in the small coupling regime the predictions of tree- level-flat inflation 
models (including many hybrid inflation models) on the size and running of the spectral 
index are quite neat and model-independent. Forthcoming WMAP analyses could 
either confirm or falsify them (e.g. if the crossing of n through n — 1 is verified). 
3 A detailed study with thresholds is performed in subsection 4.1. 
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2.2 Not-so-small-coupling regime 

If the /^-functions are positive (as usual, see the above discussion), couplings grow in the 
ultraviolet and there will be a scale where the second and third terms within the square 
brackets in (16) compete with the first one. Since they naturally have the opposite sign, 
one can expect that at sufficiently high scales (which means initial stages of inflation) 
the sign of rj, and thus n — 1, may get positive, which therefore could account for the 
preliminary indication of WMAP. Notice in particular that, since eq. (18) for rj is not 
a valid approximation anymore, the previous correlation (20) between the signs of V 
and rj does not hold, which allows rj to pass from positive to negative values without 
the presence of an extremal point of the potential V. As a matter of fact, we find 
this behaviour absolutely natural, and even unavoidable in this regime. However, to 
reproduce the WMAP preliminary indication on the running spectral index and the 
other physical requirements at the same time (in particular, a sufficient number of 
e-folds) is not a trivial matter. This will be analyzed for typical examples in sect. 3. 

The running of the spectral index can be computed, taking into account that, as 
long as f3/p <C 1, still e <C \q\ (except within the close neighborhood of the point where 
rj vanishes), so n ~ 1 + 2i], and 

d\nk- 2 d\nkd<j>~ 2MpV2e #' (25) 

where dr]/d(f> can be straightforwardly evaluated using the whole expression for r\ from 
(16), taking into account the implicit dependence on <p through the one on Q with 
the help of eq. (14). Similarly, the number of e-folds is computed using eqs. (2, 16). 
In the next section we will show explicit expressions for these quantities in particular 
examples. 

3 An explicit example 

In the previous section we claimed that, within the hybrid inflation framework, it is 
natural to expect a running spectral index, even crossing the n = 1 value, if the scale 
of the initial stages of inflation is high enough. Now we will explicitly show that this is 
the case by working out a popular example of hybrid inflation, namely the first example 
of hybrid D-term inflation, proposed by Binetruy and Dvali [15]. 
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The model is globally supersymmetric with canonical Kahler potential, and with 
superpotential 4 

W = <!>(\H + H_ - /i 2 ) , (26) 

where A and [x 2 are real constants and $, H + , H_ are chiral superfields. H± have charges 
±1 with respect to a U(l) gauge group with gauge coupling g and a Fayet-Iliopoulos 
term £ D . The associated tree-level scalar potential is given by V = V F + V D , with 

2 



V F = \\H + H_-n 2 \ + A 2 (\H_\ 2 + \H + \ 2 ) |$| 2 , 
Vd = 9 j (\H + \ 2 - \H_\ 2 + £ D ) 2 , (27) 



(the fields in these expressions are the scalar components of the chiral superfields). The 
global minimum of the potential is supersymmetric (i.e. Vf = 0, Vd = 0) and occurs 
at $ = 0, \H±\ 2 = (t£d + ^ + 4/x 4 /A 2 )/2. But for large enough |$| 

l<f| 2 >l^| 2 = ^A, A^V^+AV, (28) 

the potential has a minimum at H± = and is flat in $. This flat direction is two 
dimensional because $ is complex, $ = (<f> R + i(pj)/^/2. In the following, we examine 
the potential along the real part <fi = <pR. All our results would be the same along any 
other radial direction. It is in this flat region where hybrid inflation takes place. The 
tree-level potential during this period is simply 

v = p = ^ + \g 2 e D ■ (29) 



The one-loop correction is given by 

~ 2 2 " 

m. m. 

(30) 



~ 2 2 
~ 4 l m i 4 l m i 

In — - m 4 In — 



where m 2 (m 2 ) are the ^-dependent square masses of the fermionic (scalar) components 
of the H± superfields: 

ml = m 2 (0) = Ia 2 2 , m 2 = m 2 (0) ± A . (31) 

In practice, m\ is dominated by the m 2 (0) contribution. To see this, notice that the 
condition m 2 (0) = A (i.e. m 2 = 0) precisely corresponds to the point = C , below 

4 We follow here the formulation with a tadpole for $ in the superpotential as in rcf. [14], where 
this model was analyzed in connection with the issue of decoupling. 
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which the flatness of the potential breaks down and inflation cannot continue. Since 
the number of e-folds goes as ~ 2 , most of the e-folds have been produced at larger 
(usually much larger) values of 2 , where m 2 (0) ^> A. This is particularly true in the 
region of interest for us, i.e. the initial stages of inflation. Furthermore, as we will see 
below, inflation ordinarily ends before 2 = 2 , namely when 77 becomes 0(1), which 
normally occurs for much larger 2 . In consequence, AVi can be written as 

AVi = ^A 2 ln^M + o(-^r) , (32) 
8tt 2 Q \m\(t>)) ' 1 ; 

so V — Vo + AVi is of the generic form shown in eq. (11). 

The various /^-functions, defined as derivatives with respect to In Q, are given by 

P = ^A 2 , (33) 

A = g^ 3 > (34) 
fa D = 0, (35) 

(3 X = J_A(3A 2 -V), (36) 



16vr 2 

= i^AV- (37) 



In the leading-log approximation the radiatively-corrected potential is simply given by 
p evaluated at the scale Q = m(0) = A0/\/2; in the notation of sect. 2, 

c = X/V2 . (38) 

Concerning the end of inflation, this occurs at cn d = Max{0 c , 0^}, where C marks 
the end of tree-level flatness and is given in eq. (28), and <p, n corresponds to the point 
where r\ = 0(1). Using eq. (18) for an estimate, 

0? ^Ky (39) 

It is easy to see that C < (p v for a wide range of {g, A} so that, normally, inflation 
ends when 77 = 0(1). E.g. for g ~ A and using eqs. (3) and (7), taking into account 
that (j)(ko) in (3) and (7) should be at least equal to (f) v , one finds that this occurs for 
g, A > 9 x 10~ 2 . More realistically, when <f)(k ) ~ M p , this occurs for g, A > 10~ 3 . 

Now, a successful inflationary process, consistent with the WMAP constraints, 
should satisfy the following requirements 
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i) At least 50-60 e- folds of inflation. 

ii) A correct power spectrum of the scalar perturbations at k = 0.002 Mpc -1 , as given 

by eqs. (3) and (7). 

in) Running spectral index with n — 1 passing from positive to negative during the 
first few e-folds of inflation (after £*), with dn/d\nk = 0(— 0.05). 

As discussed in the Introduction, the experimental status of requirement in) is still 
uncertain, but we will examine here the possibilities of fulfilling it in this context. 
Likewise, condition i) can be relaxed if there are several steps of inflation, as discussed 
in footnote 1. In the latter case 7-10 e-folds (i.e. the range covered by WMAP with 
the optional complement of additional observations) might be acceptable. 

If g and A are small, (5 remains essentially constant (and positive) during inflation. 
Then the results of sect. 2 apply, in particular n < 1, with a slight running dn/dlnk < 
0, as given by eqs. (19), (22) and (23). According to that discussion, it is not possible 
in this regime to fulfill condition Hi), i.e. a running spectral index crossing through 



In consequence we will focus for the rest of this section in the case of more sizeable 
coupling constants. Let us first verify that the slow roll conditions are fulfilled. In this 
regime e and r\ are deduced from eqs. (16). Taking into account that a = 0(1) and 
that j3 is at most 0{(3), it is clear that rj = <D(M 2 f3/p(j) 2 ) and e = 0(rjP/p). Therefore 
e, \rj\ <C 1 as long as 



n — 1. 



(3 2 
— <c — — 
P Ml ' 



(40) 



which naturally leads to f3 <ti p. A simple way of guaranteeing this is to take 



A 2 9 2 e D 2 



(41) 



i.e. (3 is dominated by the Yukawa contribution: 



(3 ~ 




(42) 



8tt 2 



and p is dominated by the D-term (Fayet-Iliopoulos) contribution: 



P 




(43) 
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Imposing (41) is by no means the only way of satisfying the slow-roll conditions [e.g. 
reversing the inequalities in (41) would also work], but it is particularly simple and 
useful for a clear analytic discussion 5 . 

Using this dominance of the Yukawa coupling, we can evaluate e and rj from eqs. (16) 

2 



2 ^ \p, 



a 



02 p 



A 2 / 9 A 2 



" 1 + 2^l a+ 322^ a2 . 



(44) 



with 



As announced in the previous section, there is a scale, say Qq, at which r\ ~ 
and n — 1 passes from positive to negative. From (44), the corresponding value of A, 
A = A(Q ), turns out to be independent of the other parameters of the model, namely 

A 2 1 



— (7-^)^ 0.078. (46) 



(4tt) 2 15 

The coupling Ao is quite sizeable (A — 3.51), but well within the perturbative range, 
A 2 ,/ (47r) 2 « 1. It is perfectly possible that Q is crossed during the first few e-folds 
after t* (recall Q = \<fi/y/2), so that this hybrid inflation model (and in general any 
tree-level-flat inflation model) is able to implement the crossing of the spectral index 
through n — 1 at small k. However, to reproduce the suggested value of dn/dlnk and 
a sufficient number of e-folds at the same time [i.e. conditions i), in) above] is not that 
easy, as we show next. 

From eq. (2), the number of e-folds between Q and the end of inflation can be 
written as 



5 Actually, once global supersymmetry is promoted to supergravity, requiring /i 4 <C .g 2 £fj/2 [i.e. the 
second inequality of (41), and thus eq. (43)] is mandatory for almost any form of the Kahler potential, 
K. The reason is the well-known -q— problem, i.e. the appearance of an 0{H 2 ) mass term for the 
inflaton (and thus r\ ~ 1) if Vf ~ V. However, it is remarkable that for the supcrpotcntial of this 
model, eq. (26), and minimal Kahler potential, K = \$\ 2 , no such mass term is generated due to an 
amusing cancellation of contributions. 
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64tt 4 



~ 0.303 



Qi g 2 e D 



Ml 4 J 

qi g 2 e D 



M 2 
p 



/4 



\2/3 

A exp 



= 1.818 



16tt 2 



r l x w ( i 

, ' 1 / A end 



3 \ 2 
lQn 2 xJ 



16n 2 x 



dx 



QpH 2 
i4 



(47) 



where // = A*(Qo)- m the first line we have performed a change of variable — > Q 2 
with respect to the expression (2). In the second line we have changed variables once 
more, Q 2 — > x = 1/A 2 , neglecting the oc g 2 contributions to the RGEs of A and g 2 as 
well as the /i 4 contribution to p, according to assumption (43). Finally, in the last line 
we have taken into account that the contribution to the integral coming from small 
values of A is negligible, so we have replaced A en d —> 0, A by the value quoted in (46), 
and evaluated the integral numerically. The final result turns out to be a very good 
approximation, provided eq. (41) is fulfilled. Notice that, since the point Qo (i.e. where 
n = l) should be crossed after the first few e-folds of inflation, N® should be 0(50) 
[unless there are subsequent episodes of inflation, in which case it could be as small as 
0(1)]. 

On the other hand, the running dn/dlnk can be straightforwardly evaluated using 
(6). It gets a particularly simple form when evaluated at the scale Q , where rj = 0, 
1: 



n 



dn 



d\nk 



Qo 



-11.847 



'AT 
4vr 2 



Ml /4 



'QI 


g 2 e D 


-2 

= -0.33 


\QlH 2 ' 


-2 


Ml 


i4 . 







(48) 



where we have used (46) for the numerical estimate. Clearly, and dn/dhi k\q o 
depend on the same combination of input parameters (the one in square brackets); and 
the product (N®) 2 x (dn/d\nk\ Q() ) turns out to be independent of them 



dn 



dink 



~ 1.1 



(49) 



This constraint makes it impossible to fulfill dn/d In k\q o = 0(— 0.05) and N® ~ 0(50) 
simultaneously. However N® ~ 0(5) is perfectly possible, which means that the model 
can account for the indication of WMAP, provided there exist subsequent episodes of 
inflation. 
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logjJkCMpc 1 )] 

Figure 1: Scalar spectral index as a function of log 10 [fc(Mpc~ 1 )] for the inflaton potential of eqs. (29) 
and (30) for two different choices of parameters: one to get right N e and the other to get right dn/d\n k. 



The previous behaviour is illustrated by Fig. 1 which gives the spectral index n as 
a function of the scale wave-number k. The flat red curve [corresponding to g* = 0.02, 
A, = 4,/i 2 = (3.6xlO~ 4 Mp) 2 ,£ D = (1.2xl0^ 2 M p ) 2 , 0* = 0.28MJ is able to accomodate 
a large iV e ~ 57 but has very small running. The steep blue curve instead [with the 
same parameters except for /i 2 = (1.9 x 10~ 3 M P ) 2 and £ D = (3.6 x 10~ 2 M P ) 2 ] has a 
sizeable dn/dlnk ~ —0.075 but fails to give enough e-folds: N e ~ 7. 

Interestingly, there is no problem in reproducing the observed scalar power spectrum 
(7). Notice that P^ in eq. (3) is now 



P k 



4vr 2 1 {g\ 



2c2 \2 
Dl 



3 al\l My 



M 2 
p 



4.14 x 10 



- mi? 



~Ql 


g 2 e D 


H 4 

= 0.089^ 




Ml 


Mo . 




Mo 



(50) 



Clearly, for any given value of the quantity in the square brackets, it is possible to 
vary the parameters so that the prefactor of (50) changes in order to fit the observed 
Pk- Incidentally, the pivotal scale in the analysis of ref. [3], ko = 0.002 Mpc -1 (which 
corresponds to ~ 3 e-folds) approximately coincides with the stage at which n ~ 1 (i.e. 
when Q = Q ). This is a fortunate circumstance to perform rough estimates. 
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In summary, the simple hybrid inflation model of ref. [15], in the small coupling 
regime, produces the general results shown in subsect 2.1, with very clean predictions, 
in particular n < 1 with a mild dn/dlnk < running. In the regime of sizeable 
couplings, the model can account for the scalar power spectrum at the pivotal scale 
and a running spectral index crossing n — 1 at a suitable scale with the suggested 
slope, dn/d\nk = 0{— 0.05). However, the number of e-folds produced is rather small, 
so this possibility is not viable unless there are subsequent episodes of inflation. 

The previous analysis can be repeated for an arbitrary number of pairs of H± fields, 
using the formulae given in Appendix A. The results do not change much, except for 
the fact that the value of A scales as 1/y/N. Also, the right hand side of eq. (49) 
depends slightly on the value of N, being always ;$ 1.1. 

Of course, one cannot exclude that other tree-level-flat models are able to reproduce 
such running spectral index without the above limitation for the number of e-folds. To 
achieve this one needs sizeable r\ at high scale, so that the spectral index runs apreciably 
in the the first stages of inflation, and small e at lower scales, so that a sufficient number 
of e-folds is produced. This implies that (3 (and thus the relevant coupling constants) 
must evolve from sizeable values at high scale to substantially smaller values at lower 
scales. The faster such evolution, the easier to accomplish both requirements. Since this 
evolution is a result of the RG running, one can imagine scenarios where it is speeded 
up. For example, if some threshold of new physics is crossed in the inflationary course, 
the /3-functions may suffer quick changes, as desired. However the analysis has to be 
done carefully to eliminate spurious effects, as we shall see shortly. On the other hand, 
the existence of new physics can have other (even more) interesting implications. 

In the next section we analyze carefully the various effects of thresholds of new 
physics in the inflationary process. 

4 Thresholds of new physics 

The Lagrangian responsible for the inflationary process can well be an effective theory 
of a more fundamental theory at a higher scale. Actually, since during inflation the 
characteristic scale is changing (following the value of the inflaton field), it may even 
happen that some threshold of new physics, say A, is crossed in that process. In either 
case one can expect effects from the new physics in the details of inflation. 
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For the kind of inflationary models we are considering in this paper, namely tree- 
level-flat models, there can be two main types of new physics: 1) new physics that 
does not spoil the tree-level flatness of the potential and 2) new physics that does 
spoil it. In case 1), the new physics can only affect the form of the initial Lagrangian 
through radiative corrections. If A is above the scales of the whole inflationary process, 
these effects are suppressed by the ratio of the scales, ~ (0/A) n , besides the typical 
loop suppression factors. Thus we expect them to be much less important than the 
radiative effects calculated within the effective theory considered in previous sections. 
On the other hand, if A is crossed in the inflationary course, the effective potential and 
the various RGEs involved will get modified, which can have interesting implications, 
as mentioned at the end of sect. 3. In contrast, for case 2) the new physics can show 
up, even below A, through non-renormalizable operators in the tree-level potential, 
suppressed by ~ (0/A) n . Actually, in case 2) it is not realistic to suppose that A is 
crossed during inflation, since above A the new physics will usually spoil the slow-roll 
conditions. Next we will explore possibilities 1) and 2) in turn. 

4.1 Thresholds crossed in the inflationary process 

Let us suppose that the inflaton potential is flat at tree-level and the 1-loop radiative 
correction, AVi, contains contributions from states with different masses. This gener- 
alizes the situation discussed in sect. 2. To be concrete, consider two types of particles 
with masses m((f>) and M((f>) satisfying m((f>) <C M((f>). The 1-loop effective potential 
reads 

W)=M Q) +A l„I^ + ( A - fl ), n £ffi. (51) 

The tree-level potential, p, depends implicitly on Q, as indicated, with RGE cor- 
responding to the full theory, dp/d\nQ = (3 h . Notice that f3i is the corresponding 
(3— function when only the light particles are present. Generically, m 2 (0) = m 2 + c 2 2 , 
M 2 (4>) = M 2 + C 2 (f) 2 , where m and M do not depend on 0, and c and C are coupling 
constants. 

As mentioned above, we focus on the interesting situation in which the scale of 
inflation crosses the threshold M where the heavier particles decouple. In that case 
one needs a reliable prescription to evaluate V(<f>) above and below the threshold. In 
particular one should take care of potentially large logarithms in the potential (51) 



16 



and higher order corrections to it. In the case of just one type of mass, say m(<f>), this 
can be achieved by choosing Q = m((f>), so that the unique log in the radiative cor- 
rection cancels and V(<p) = p[Q(<f))]. This corresponds to summing up the leading-log 
contributions in V to all loops. This simple procedure is not possible here because 
of the two different logs in (51). Then, a convenient approach is to perform the inte- 
gration of the RGEs in an effective theory framework. To this end, we consider two 
different regions for the potential. In the high energy region, when {Q,m((f))} > M, 
the potential, Vhi g h, is as written in (51), with RGEs including the virtual effects of 
the heavy states. However, in the low energy region, defined by {Q,m((f))} < M, the 
potential, Vi ow , is defined by dropping the contribution from the heavy multiplets to 
the potential and to all RGEs. For this procedure to be consistent one has to match 
the high- and low-energy potentials, i.e. one has to add a piece 5thVi ow (4>) to V\ ow to 
guarantee Vi ow = Vhigh at Q = M. In this way one gets 



where the subindex M means that the quantity within the parenthesis is evaluated at 
Q = M. The additional piece S t hVi ow ((f)) can be expanded in powers of <j)/M if desired: 
one would get <fi operators suppressed by inverse powers of M, which is the reason one 
can discard such high energy remnants when interested in low energy physics. However 
we keep here the whole expression as we are also interested in the potential for values 
of 4> in the neighbourhood of the threshold. Actually a correct treatment of that region 
is crucial for a reliable analysis. Notice for instance that since e ~ 1/0 most of the 
e-folds in the energy region below Q = M are produced precisely near that threshold. 
The inclusion of the 5 t hVi ow (<f)) contribution usually increases e, therefore decreasing 
(importantly) the number of low-energy e-folds, and should not be neglected. 

The next step is to make a judicious choice of the renormalization scale Q. The 
simplest option is to take Q((f>) = m{4>), which works well in both energy regimes 
(meaning that the explicit logarithms are never large). After doing that we get the 
following expressions for the potential in the two regions: 




(52) 




(53) 



V h igh(4>) 



p[Q(0)] + In 



M(0) 



m{4>) ' 
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v low (<f>) = pIQ(<P)} + (p h - A) \ M ln 



M 



(54) 



As a matter of fact, the previous approximation may not be accurate enough in 
some cases, e.g. when the coupling constants involved are not so small. In particular, 
although V in (54) is continuous across the threshold by construction, the derivatives 
V, V" (and thus e and if) are not. This is an spurious effect which is conveniently 
smeared out when the approximation is improved. This can be done by including 
higher-loop corrections to the effective potential. Indeed, without any further calcula- 
tions, one can obtain the leading-log corrections at arbitrary loop order by using the 
Q-invariance of the potential. In particular, at 2-loop leading-log (2LL) the potential 
reads 



V 2LL ((j)) 



p(Q) + Aln^ + (^-A)ln M(0) 



Q 



Q 



In 



m(4>) 



+ -[(3 h -2 (3 l +A 



In 



M{4>) 

Q 



+ Pi -A 



(55) 



Q Q 

where the circle (dot) indicates d/dlnQ in the high (low) energy region. Redoing 
the decoupling and matching procedure at the same scale Q = M, the corresponding 
expressions of V 2LL in the high and low energy regions read 



p[W)] + (/^-A)m 



m(0) 



+ I 1/3, -2 A +A 



ln 



M(0) 



y2LL 
low 



(0) 



p[Q{<t>)\ + W h - A) U ^ ^ +\[P h -ih +A 

ln m(0) ^ M(0) 









M 



+ [Pi -Pi 



(56) 



m M M 

The V, V" derivatives are straightforward to evaluate by taking into account the 
0— dependence of the scale, i.e. Q 2 {4>) = m 2 + c 2 4> 2 . An alternative procedure, starting 
with the 1-loop initial potential (51), to incorporate the leading- log corrections at 
arbitrary order is to perform the decoupling at the 0— dependent scale M(0) (which 
leaves no threshold corrections) and then evaluate the potential at the m(4>) scale. In 
this way one obtains expressions for the potential in the high and low energy regions 
which coincide with those of the previous procedure. In particular, collecting the 2-loop 
leading- log contributions one exactly recovers the result (56). 
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To illustrate this general scheme of constructing the effective potential across thresh- 
olds of new particles, let us consider again the simple hybrid inflation model of sect. 3, 
but now with iV additional pairs of {H + ,H^} fields with invariant mass M. The 
superpotential of the full (high-energy) theory reads 

/N+l \ N+l 

W = $ {J2 XH l H - - A* 2 + MH+H" , (57) 

\a=l / a=2 

where, for simplicity, we have taken equal Yukawa couplings for all the {H + , HJ\ pairs. 
For large enough |$| the tree-level potential is flat, having the same form as in eq. (29) 

V = p = p 2 + l -g 2 il . (58) 

The 1-loop effective potential has the form (51) (up to suppressed corrections) with 

m\<P) = \\ 2 <P\ M 2 (0) = M 2 + iA 2 2 , (59) 



where = y^Re $ and 



A = ^A 2 , ^ = ^(iV + l)A 2 , (60) 



with A = sjg^l + A 2 /i 4 - 

Using the above RGEs, we can also include the two-loop leading log corrections, 
simply making use of the general expressions (51) and (55): 

V(<f>) = P (Q) + -^A 2 (L m + NL M ) 

+ (8^)2 { A V [ 2L ™ + 2NL ML m + N(N + 1)L 2 M ] - 2\ 2 g 2 p\L 2 m + NL 2 M ) 
+ g 6 &(2L 2 m + NL m L M + N 2 L 2 M )} , (61) 

where we have used the short-hand notation L m = \n[m(<fi) / Q] and L M = ln[M(0)/Q]. 
We have checked that this potential agrees with the one obtained by using the ex- 
pressions for the two-loop potential of generic supersymetric models given in ref. [16] 
(which can also be used to add subleading logs and finite corrections). 

The expressions for the potential above and below the threshold at M can be easily 
obtained following the effective approach described above. More precisely, one can 
write the potential in the high- and low-energy regions up to 1-loop or 2-loop leading 

o 

log order using eqs. (54) and (56) respectively. The $i t h and @ l h derivatives for the 
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latter can be easily extracted from (61) [by comparison with (55)], or calculated using 
the /9-functions for the various parameters of the model involved in the expressions 
(60), which are given by eqs. (34-37) for the low-energy region and by eqs. (A.4-A.7) 
for the high-energy one. We do not present explicitly the corresponding formulas. 

The results are somewhat better than in the case without thresholds analyzed in 
sect. 3. Nevertheless the conclusions are basically the same: The number of e- folds and 
the slope of the spectral index are correlated by an equation similar to (49), so that 
it is not possible to get dn/d\nk = 0(— 0.05) and N e = 50 — 60 at the same time. It 
is worth remarking that this happens in spite of the fact that the f3 functions undergo 
quick changes in the threshold region, as desired. However, the corrected potential 
[e.g. as given in eq. (56)] softens these effects and the net impact in N e and dn/dlnh 
gets much reduced. 

4.2 Thresholds of new physics above the inflationary scale 

We consider now an even simpler modification of the inflationary potential. It assumes 
a scale of new physics, M, higher than the scales relevant to inflation. Nevertheless, 
inflation is affected through the presence of non-renormalizable operators of the inflaton 
field that lift the flat direction along which the inflaton rolls. It could appear naively 
that such operators should have a negligible impact on inflation whenever they can be 
reliably taken into account (i.e. <p 2 /M 2 <C 1) but we show in this section that this is 
not the case. In fact, through such effects we are able to give a simple and complete 
potential that satisfies all three goals listed in section 3. The analysis can be 

carried out in a very model-independent way. 

Let us start directly with an inflaton potential that reads 

VW=p + ^Hf. + ^. (62) 

The first two terms just correspond to the generic one-loop potential we have discussed 
in previous sections. In the small-coupling regime, as in subsect. 2.1, we take (3 as a 
constant. The last term in (62) is a non-renormalizable operator (NRO) left in the 
low-energy theory after integrating out some unspecified physics at the high scale M. 
This scale absorbs any possible coupling in front of the operator. Of course V(<f>) may 
contain other NROs of different order. Here we assume that the one shown in eq. (62) 
is the lowest order one; NROs of higher order will have a negligible impact compared 
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to it (provided they are suppressed by the same mass scale M). The sign and power we 
have assumed for this NRO are convenient to guarantee the stability of the potential. 
Notice also that an even power for this operator is what one expects generically in 
super symmetric theories. In fact, the question of whether such potential can have a 
super symmetric origin is very interesting and we leave its discussion for the end of this 
section. 

Now we focus on studying whether a potential of the form (62) can reproduce the 
WMAP indications for a running n, and for what values of the potential parameters 
can it do so. By trivial inspection of the derivatives of V(4>) with respect to 0: 



we realize that the NRO can have a significant impact on inflation when the small 
number (4>/M) 2N is comparable in size to /9/0 4 (which is also quite small). In other 
words, even for <fi 2 <C M 2 , the NRO can have an important effect because it gives a 
correction to a potential which is almost flat. 

It is also immediate to realize from (63) that, for sufficiently large 0, the higher 
derivatives V", V" (and thus n, £) can receive a large contribution from the NRO while 
the contribution to V (and thus e) is much less significant, thanks to the additional 
(2N + 3) and (2N + 3)(2iV + 2) factors in V", V" . In consequence, as <fi rolls down its 
potential the effect of the NRO quickly dies away and inflation proceeds as discussed 
in section 2 but, in the early stages of inflation the running of n can get important 
modifications from the NRO corrections, even though the values of V and (to some 
extent) V are scarcely modified. Let us examine this in more quantitative terms. 

The slow-roll parameters (1) can be readily found as a function of from the 
derivatives of the potential listed in (63). The e parameter is much smaller than rj and 
£, so n ~ 1 + 2i] and dn/d\nk ~ — 2£. In this way we obtain 



V'"{<t>) 



v'\<t>) 



v'(<f>) 




M 2N ' 



(63) 



dn 



1 



[l + (2iV + 3)(iV + l)r(0)] [l + r(0)] , 



(64) 



dink 



N 2 M 



where 



r{<j>) = 2{N + 2) 



4 <p 2N 

~Jm 2N ' 



(65) 
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Figure 2: Running of the spectral index, given by — dn/d In k\*, as a function of iV for several 
values of N°/N e . Red solid (open blue) dots have N e = 50 (60) and iV e ° = {47.5,45,42.5,40} 
({57.5, 55, 52.5, 50}), with higher values of N® giving smaller running. 



and 



N e (<f>) 



2/3M? 



(66) 



which is approximately the number of e-folds from <ft till the end of inflation [obtained 
by neglecting the NRO contribution to V'((p) in eq. (2)]. Applying eq. (64) to the 
starting point 0*, it looks possible to get N% (dn/d In ~ —125, as suggested by 
WMAP, if N and/or r(<p*) are not small. We analyze next whether this is feasible. 

Let us call the value of the inflaton at that particular point with n — 1. At 0o we 
should have V'(0o) — 0, which implies r(0 o ) — 1/{2N + 3). Using the approximation 
of eq. (66) it is simple to get 

1 {N\ N+2 
r(<t>.) ~ r \U\ , (67) 



2N + 3 VA'e / 

where N e = N e (<fi*) is the total number of e-folds, and N® = iV e (0 o ); ^e. n — 1 
after the first N e — e-folds. Using Eqs. (64-67) we can find out what values of 
iV are required to get enough running for n. This is shown in Figure 2, which gives 
dn/dln as a function of N for different choices of N®/N e [corresponding to &(0o) 
not far from 10~ 2 Mpc -1 ]. From this plot we see that large values of N are required to 
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Figure 3: Inflaton effective potential (normalized to p) with a NRO as in eq. (62) for N = 9 and 
parameters as given in the text. The star marks and the circle, (f>o- 



get a significant running of n, with smaller N®/N e being preferred. As an example, let 
us take N = 9 and N°/N e = 42.5/50, which gives dn/d\nk\* ~ -0.03. 

For a fixed value of the scale of new physics, M, the quantities j3 and p are de- 
termined by the constraints of Pf. and N e : one finds p ~ (10 _3 y / MM p ) 4 and (3 ~ 
(10 _4 M) 4 . Choosing for instance M ~ 0.95M P one further gets dn/dlnk\* ~ —0.03 for 
(f)* ~ 0.15M p while 0o — 0.142M p ; the numerical values for the number of e-folds are 
N e ~ 48.8 and A^" e ° ~ 42.4; and one gets P k ~ (2.95 x 10" 9 ) x (0.8). Note that (f)jM is 
sufficiently small for the effective theory with the NRO to be trustable. Figure 3 shows 
the effective potential as a function of <p/M p , with 0* indicated by a star and 0o by 
a circle. Notice how 0* is below the range where the NRO starts to be important for 
V'((j)) [but not for V" (</))]. The upper plot in Figure 4 shows the slow-roll parameters 
as a function of <fi/M p . Inflation does not continue below <fi ~ 0.02M p . Finally, the 
lower plot of Figure 4 gives the scalar spectral index as a function of log 10 fc(Mpc _1 ). 

The running of the spectral index with k, as shown in the picture just mentioned, 
can be well approximated by 

1 1 

n ~ n * + T e ~ N e - hpg 
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/ 1 k \ JV+1 

The first line is the same as eq. (24), which is valid for (weak coupling) potentials with 
a radiative dependence on the logarithm of the inflaton field. The effect of the NRO is 
explicitly seen as the additive term in the second line. 

As already discussed in connection with Figure 2, stronger running of the spectral 
index can be achieved if N is larger. This leads us to the next important question: how 
reasonable is it to expect a non-renormalizable operator with e.g. N = 9 or even larger? 
Recall that, for the previous mechanism to work, such operator should be the leading 
one. This is not a serious problem because supersymmetric flat directions as the one 
we are using for inflation can be protected against lifting by additional symmetries. 
In fact it is common that some of the flat directions of supersymmetric models are 
only lifted by NROs at very high order. This is well known, for instance in the MSSM 
[17]. This is also very common in the context of D=4 string theories, where string 
selection rules forbid many operators in the superpotential that would be allowed just 
by gauge invariance. For example, in the popular Z 3 — orbifold compactifications of the 
heterotic string, non-renormalizable couplings involving twisted matter fields (which 
would be the relevant case if the inflaton is a twisted field) have the structure $ 3+9ri ; 
with n = 1,2,... [18]. 

Concerning the supersymmetric realization of a potential like that in (62) the first 
guess would be to use a superpotential of the form 

1 $ N 

W = <K*H + H—f) + WT jf> wf , (69) 

which is simply the standard superpotential we used in section 3 supplemented by 
a non-renormalizable term for the inflaton. From (69), one obtains at tree level the 
inflaton potential 

(70) 

where Vd is the Fayet-Iliopoulos contribution. This potential has a term of order 
$ 4 $ 2Ar /M 2iV like the one we are after, but also a term of lower order (-(j?<fr 2 <& N /M N + 
h.c). For the higher order term to dominate one would need |/x 2 | <§C <§> 2 <& N / M N which 
can in principle be arranged. In any case, the presence of a non-zero \x poses a different 
problem. In view of the form of the potential (70) it is clear that there will be a series 



dn 



(N + l) \d\nk 




V = V D + 



N 



M N 



25 



of minima with $ = $ m exp[i7r/(2 + iV)] where $£f +2) = \n 2 \M N . The flat potential for 
$ is therefore lifted by the NRO along some directions in the complex plane of $ while 
along other directions the potential develops minima and this can change qualitatively 
the evolution of the inflaton. It is interesting that the usual logarithmic one-loop 
corrections that cause the inflaton to roll in the standard scenario can also cure the 
previous problem in the present case so that the minima are no longer dangerous. In 
that case one can have the inflaton running in a stable trajectory along which only its 
real part is non-zero. Nevertheless, the simplest cure of the previous problems is to 
choose /i — which, incidentally, is the choice of the original formulation of the D-term 
inflation model [15]. In that case one simply gets a potential term of order 4 (/> 2Ar /M 2Ar 
as in eq. (62). Notice how, thanks to supersymmetry, the order of the NRO in the 
superpotential (69) has almost doubled in the potential (70). 
An appealing alternative is to start with 

W = A<M+7/_ + Im*' 2 + jp^f^ , (71) 

where we have introduced an extra field $' with mass m < M (so that $' really belongs 
in the effective theory below M). When one is along the flat direction for the inflaton 
field $ the coupling to $' generates a tadpole for it so that also $' develops a VEV. 
As we are only interested in the evolution of the inflaton field, we can simply eliminate 
$' by solving its equation of motion in terms of <£> and substituting in the two-field 
potential V(<&, to get the final potential for $: V(&) = l / ($, $'($)). In this way 
one gets the following inflaton potential (along = a/2 Re <3> and absorbing factors of 
y/2 in M) 

where again Vb is the Fayet-Iliopoulos contribution. The advantage of this option is 
that the NROs are naturally of higher order. For instance, to get N = 9 one needs a 
rather modest P = 4 in the superpotential (71). 

5 Conclusions 

In this paper we have assumed that the WMAP indications on the running of the 
scalar spectral index are real and then explored the intriguing possibility that this 
runnig is caused by physics at very high energy scales affecting the inflaton potential. 
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It is tantalizing to imagine that we could be seeing in the sky the imprints of such 
physics. In more concrete terms we aim for well motivated inflaton potentials (from 
the point of view of particle physics) that can reproduce the suggested slope of n with 
scale, dn/dlnk ~ 0(— 0.05), giving at the same time the right amount of e-folds and 
the right amplitude of the power spectrum of scalar fluctuations. 

This is not an easy task and has atracted some effort recently. We have focused on 
a large class of inflationary models which we find particularly atractive: those with flat 
tree-level potential lifted by radiative corrections. These radiative corrections directly 
control the slow roll of the inflaton and, in turn, the running of the spectral index n. 
These models include some typical hybrid inflation models as a particularly interesting 
subclass. In these models there is a one-to-one correspondence between the value of the 
inflaton potential, the related distance scale k and the associated renormalization scale, 
so that the running of the inflaton down its potential scans different high energy scales 
and cosmological distances. In this language, WMAP indicates that the spectrum of 
scalar perturbations is blue (n > 1) at the largest cosmological scales (k ~ 10~ 4 Mpc _1 ) 
and turns to red (n < 1) below some scale k ~ 10~ 2 Mpc -1 . In terms of energy scales, 
at the start of inflation the inflaton had a very large value (corresponding to very 
high energy scales) and its slow-roll produced n > 1 and a negative slope dn/d\nk. 
After a few e-folds this negative slope made n cross 1 at some particular value 0o of 
the inflaton field, corresponding to some particular energy scale Q Q . Later on, n gets 
stabilized around 0.94 — 0.99 and the additional number of e-folds accumulates till 
A^g ~ 50 — 60, when inflation ends. 

In our search for suitable potentials in this type of models we find two different 
regimes depending on how sizeable the coupling constants of the model parameters 
are. In the small-coupling regime, these models make sharp predictions for the size of 
the spectral index n and for its running with scale, namely 



From this equation, it is clear that the preliminary WMAP indication of n crossing 
n = 1 cannot be reproduced. 

On the other hand, at larger coupling we find that a much stronger running of n is 
absolutely natural and even unavoidable. However, reproducing the WMAP indication 
for dn/d\nk and a sufficient number of e-folds (N e ~ 50 — 60) simultaneously is very 




(73) 
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difficult. We have illustrated this in a popular D-term hybrid inflation model for which 
we have obtained a 'no-go' constraint of the form (49) 

dn 



(N ) 2 

1 e) dink 



~ -1.1 , (74) 



where, as mentioned above, Q is the scale at which n crosses 1, and N® is the number of 
e-folds from that scale till the end of inflation. Now, if one requires dn/dlnk ~ —0.05, 
then N® <C 50. Therefore, unless a subsequent stage of inflation completes the total 
number of e-folds, these models are in trouble by themselves to reproduce the behaviour 
of n mentioned above. 

As the details of inflation in this region of scales are sensitive to possible physics 
thresholds, we considered the natural possibility of the inflaton coupling to heavy new 
fields with a mass M that falls in the energy range crossed by the inflaton field in the 
first few e-folds. Somewhat to our surprise we found that, in the presence of sharp 
thresholds, the results of the previous case are not changed substantially and still one 
finds that the requirement of a significant running of n leads to a small number of 
e-folds being produced. 

A different way in which high energy physics can leave an imprint in the predictions 
of inflation is through non-renormalizable operators that modify the inflaton potential 
in the inflationary energy range . In this case, the new energy threshold is not crossed 
during inflation (this is necessary for a reliable treatment of new physics effects in terms 
of effective operators in the first place!). The key point here is that even operators that 
are suppressed by a high energy mass scale can become important in the environment 
of a very flat potential. In this context we are able to write down very simple and well 
motivated inflationary potentials that can give a running spectral index in agreement 
with WMAP and a sufficient number of e-folds for reasonable choices of the parameters. 
The succesful type of potential combines three ingredients: tree-level flatness, smooth 
dependence on the inflaton through radiative logarithms and a non-renormalizable 
term: 

V W =p + /3 .„^ + ^. (75) 

We have further motivated this type of potentials by finding supersymmetric realiza- 
tions through superpotentials that include non-renormalizable operators (NROs) of 
modest order. 
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Finally, we summarize the functional dependence of n(k) expected for flat tree-level 
potentials. In the small-coupling regime, and if no mass thresholds are crossed during 
inflation, this is simply given by 

B( * ) = 1 - *.-M*A.r (76) 

where fc* is the initial distance scale ~ 1CT 4 Mpc -1 . The only parameter in (76), 
N e , corresponds typically to the total number of e-folds, i.e. N e ~ 50 — 60, although 
it could be either smaller [if there are subsequent episodes of inflation], or larger -or 
even negative- [if the end of inflation is not marked by r] = 0(1), but by the violation 
of some hybrid-inflation condition for the tree-level flatness]. This is explained in 
subsect. 2.1. The case with thresholds and/or not-so-small couplings are addressed in 
subsects. 2.1, 2.2, 4.1. 

On the other hand, when there are effects of high-energy (new) physics through 
non-renormalizable operators (NROs), n(k) is of the form 

1 1 



n(k) = n* + 



N e N e -ln(k/K) 



N e ( dn 



(N + 1) \d\nk 



+ N? 



■lln*Y +1 
N P n kJ 



(77) 



where * denotes quantities evaluated at the initial scale k*, and N is the exponent of 
the dominant NRO, as expressed in eq. (75). The quantity N e has the same meaning 
as in eq. (76). The other parameters in eq. (77) are n*, dn/dhik\*, which can be 
expressed in terms of the parameters of the inflaton potential (75), as explained in 
subsect. 4.2. Finally, if the effect of the NROs is negligible [e.g. if M — > oo in (75), 
which corresponds to n* — > 1 — 1/N e , dn/d\ak\* = —1/N^ in (77)], then eq. (76) is 
recovered. For a more detailed discussion, see subsect. 4.2. 

It would be a very nice test for the wide and well-motivated class of inflation models 
analyzed in this paper, if future analyses and measurements of the spectral index n(k) 
could be adjusted by one of the previous expressions, which in turn would provide 
precious information about the inflaton potential and on physics at very high energy- 
scales. 
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A. The simple hybrid- inflation model with N+l flavours 



We consider here the same hybrid- inflation model of sect. 3, with N + 1 pairs of H± 
fields, instead of the unique pair (i.e. N = 0) of the original model. The results have 
interest on their own (e.g. they illustrate how the values of the coupling constants and 
the speed of running are affected by the number of flavours), but they are also useful 
to illustrate the threshold-crossing procedure discussed in sect. 4.1. 
The superpotential reads now 



JV+l 



(A.I) 



a=l 



where A a are the N + l Yukawa couplings. As for the A^ = case, H± have charges 
±1 with respect to the U(l) gauge group with gauge coupling g and Fayet-Iliopoulos 
term £ D . The tree-level scalar potential is given by V — V F + V D , with 

2 



Vp = 



V, 



D 



E(V: 



9_ 
2 



+ £A^|tf_| 2 + |// + | 2 )|<t>| 2 , 



(A.2) 



Again, the global minimum of the potential is supersymmetric (Vp = 0, Vd = 0), but 
for large enough |$| the potential has a minimum at H± = and is flat in 0. The 
tree-level potential in this region is 

1 



The various /3-functions, defined as derivatives with respect to In Q, are given by 

1 



Pin 



8tt 2 
0, 

1 



(N+l)g* , 



16tt 2 



3A 2 + E A '- V 
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2 E ^ , 
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Af+1 

(N + l)g A & + E W 
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167T 
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87T 



(A.3) 

(A.4) 
(A.5) 

(A.6) 

(A.7) 
(A.8) 



For simplicity we can take the same Yukawa coupling A a = A for all pairs, 
which is a stable condition under RG evolution. 
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